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BINARY-BINARY INTERACTIONS 
AND THE FORMATION OF THE 
PSR B1620-26 TRIPLE SYSTEM IN M4 

Frederic A. Rasio 1 ' 2 , Steve McMillan 3 , and Piet Hut 1 
ABSTRACT 

The hierarchical triple system containing the millisecond pulsar PSR B1620— 26 
in M4 is the first triple star system ever detected in a globular cluster. Such 
systems should form in globular clusters as a result of dynamical interactions 
between binaries. We propose that the triple system containing PSR B1620— 26 
formed through an exchange interaction between a wide primordial binary and a 
pre-existing binary millisecond pulsar. This scenario would have the advantage 
of reconciling the ~ 10 9 yr timing age of the pulsar with the much shorter lifetime 
of the triple system in the core of M4. 

Subject headings: celestial mechanics, stellar dynamics — globular clusters: gen- 
eral — globular clusters: individual (M4) — pulsars: individual (PSR B1620-26) 
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1. Introduction 

The millisecond pulsar PSR B1620— 26 in the globular cluster M4 has a low-mass com- 
panion (most likely a white dwarf, of O.3M for a pulsar mass m p = 1.35 M & ) 
in a nearly circular orbit of period Pi = 0.524 yr (Lyne et al. 1988; McKenna & Lyne 1988), 
corresponding to a separation a\ % 0.8 AU. The very large value of the pulse period second 
derivative P = -2.3 x lO" 2 ^" 1 (Backer 1993), indicates the presence of a second, more 
distant orbital companion (Backer 1993; Backer, Foster, & Sallmen 1993; Thorsett, Arzou- 
manian, & Taylor 1993). The current timing data are consistent with a second companion 
mass in the range m 2 ~ 10~ 3 — 1 M@, with a corresponding semi-major axis a 2 ~ 10 — 100 AU 
for the outer orbit (Michel 1994a, b). Recently, a candidate optical counterpart consistent 
with being a main-sequence (hereafter MS) star of IXlclSS 777-2 ~ 

0.5M 9 has been identified 

by Bailyn et al. (1994). In this paper we assume that this optical identification is correct 
(and therefore we do not consider the possibility that the second companion might be a 
Jupiter-like planet; see Sigurdsson 1993), and we explore the consequences theoretically. 
Our results would not be changed significantly if the second companion were a MS star of 
somewhat lower mass. A MS star more massive than 0.5M 9 is strictly ruled out by the 
optical observations. 

The binary pulsar's orbital eccentricity, e\ = 0.025, although very small, is several or- 
ders of magnitude larger than observed in most other low-mass binary millisecond pulsars 
(Thorsett et al. 1993). If this anomalous eccentricity has been induced by secular pertur- 
bations in the three-body system, then the second companion mass must be ~ 0.1 — 1 M & 
(Rasio 1994a, b), consistent with the optical identification of Bailyn et al. (1994). Alterna- 
tively, the eccentricity might have been produced by the dynamical interaction that formed 
the system. We discuss this possibility briefly in §2. 

Hierarchical triple systems are expected to be produced quite easily in globular clusters 
through dynamical interactions involving primordial binaries (Mikkola 1984; Hut 1992). 
There is growing evidence that the primordial binary fraction in globular clusters may be 
as large as 30%, with 10% generally considered "typical" (Hut et al. 1992b). When binaries 
dominate the core population (because of mass segregation), a few percent of them are 
expected to be in triples (McMillan, Hut, & Makino 1991). Here we consider in detail 
interactions between a wide primordial MS binary and a binary millisecond pulsar (hereafter 
BMP). In a significant fraction of these interactions, one of the two MS stars is ejected 
while the other remains bound in a wide orbit around the BMP. We note that interactions 
between a binary and a single star cannot produce a stable hierarchical triple in the absence 
of dissipation (Chazy 1929; but see Bailyn 1989). 
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Our motivation for invoking a pre-existing BMP is the following. If the second pulsar 
companion is indeed a MS star of O.5M , then the eccentricity of the outer orbit 

must be large (e 2 > 0.5) and its semi-major axis cannot be much smaller than ~ 50-100 AU, 
to be consistent with the pulsar timing data (Michel 1994a, b; Rasio 1994a, b). Such a wide 
orbit is easily disrupted in the dense core of M4 (at present PSR B1620— 26 is most likely 
just outside the core but the density there is not significantly lower): its lifetime is then 
~ 10 7 — 10 8 yr (see §4). This is much shorter than the age of the BMP, which is most likely 
^ 10 9 yr (see Thorsett et al. 1993). Thus the millisecond pulsar could not have formed inside 
the triple system. 

In §2 and §3, we present the results of detailed numerical calculations of interactions 
between a BMP and a wide primordial MS binary. These calculations were done with the 
automated software package STARLAB (Hut 1994; Hut & McMillan 1994). The long-term, 
high accuracy triple integrations described in this paper are made possible by the use of a new 
time-symmetrized integration algorithm, whose remarkable stability properties are described 
in more detail by Hut, Makino, & McMillan (1994). In §4 we discuss the formation rates of 
triple systems like PSR B1620-26 . 

2. Eccentricity Perturbations 

Since the eccentricity of the inner BMP in the presently observed system is very small, 
the BMP cannot have been perturbed very much during the exchange interaction 4 . This 
places a limit on how close either of the two MS stars can have approached the BMP. In our 
calculations of exchange interactions (§3), the BMP is treated as a single mass and given 
an effective radius R e such that close passages with periastron separation r p > R e cannot 
produce an eccentricity larger than observed today. We first determine an appropriate value 
for R e . 

Fig. 1 shows the induced eccentricity Aei as a function of pericenter distance. We assume 
that the BMP had a perfectly circular orbit before the encounter and we place the perturber 
on a slightly hyperbolic orbit with relative velocity at infinity ~ 8kms _1 , comparable to 
the 3-D velocity dispersion in the cluster 5 . We see that the constraint Aei < 0.025 (the value 



4 Tidal circularization is completely negligible given the large binary separation oti 0.8 AU. Even if the 
pulsar's inner companion were a fully convective low-mass MS star, the eccentricity damping time would be 
~ 10 18 yr (using eq. [6.2] of Zahn 1977). In reality, the inner companion is most likely a white dwarf and the 
effect is even smaller. 

5 Because we are considering marginally soft primordial binaries, this is also the typical orbital velocity of 
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observed today) implies that no star of mass 0.5 M & can approach within R e % 3«i ~ 2.4 AU 
of the BMP. This decreases to R e % 2.5«i for a perturber of mass O.2M and increases to 
R e ~ 3.5ai for O.8M . In addition, for random orientations and phases, the distribution 
of Aei values around the average is rather wide (and markedly non-Gaussian, with a long 
tail extending to small values of Aei; see Fig. 1, insert). Thus one cannot strictly define an 
effective radius R e that applies to all situations. However, given the enormous computational 
advantage of treating the BMP as a single mass, we nevertheless adopt this prescription as 
an approximation. 

It is possible that a close encounter during the exchange interaction might have induced 
the eccentricity of 0.025 observed today. However, we find this possibility less attractive for 
the following reasons. First, because the dependence of Aei on r p is very steep (close to 
the Aei ex exp(— r^ 2 ) relation expected for small deviations from adiabatic invariance; see 
Phinney 1992, Heggie & Rasio 1994), one would have to tune the parameters quite carefully 
to obtain an eccentricity near the desired value. In addition, as shown by Rasio (1994a, b), a 
second companion of mass m 2 ~ 0.5 M & can induce very naturally this measured eccentricity 
through secular orbital perturbations, so there is no real need for any other explanation. 

3. Exchange Interactions 

We now consider interactions between the BMP and a MS binary, and we determine the 
cross section for triple formation and the typical characteristics of the triples formed through 
this mechanism. We represent the BMP as a single mass with effective radius R e = 3«i. 
One of the two MS stars must be the present outer member of the triple, of O.5M . 
The other can have any mass but is most likely of mass <^ O.5M since the least massive 
component is most likely to be ejected during the interaction. For definiteness we have 
considered the three values 0.2, 0.5, O.8M . 

Fig. 2 gives the total interaction cross section E as a function of clms, the semi-major axis 
of the incoming MS binary. Each cross section was determined by simulating numerically 
~ 5 X 10 4 interactions. The statistical error on E is always less than 5%. Energy equipartition 
and a Maxwellian distribution for the relative velocity v r have been assumed. The average 
(Etv) shown in Fig. 2 is over this Maxwellian distribution. We adopt the value a = 5kms _1 



a MS star during the exchange interaction. We did not attempt to perform an average over relative velocities 
since the distribution of orbital velocities during the interaction is not known, and, in particular, need not 
be Maxwellian. The results are not sensitive to the adopted value for Voo, since the trajectories are never far 
from parabolic. 
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for the 1-D velocity dispersion of O.8M objects. For interactions between objects of mass 
m a and nib and 1-D velocity dispersion a a and o^, the relative (3-D) velocity dispersion 
is then given by a 2 = 3(cr 2 + a%) with m a a 2 = m^al = 0.8 M & a 2 . The eccentricity of 
the primordial binaries is assumed to follow a thermal distribution, with probability density 
p(e) = 2e (Heggie 1975). We see that in the region of interest we have (Y, tot v r ) ~ AR e aMS&r, 
with the constant A % 200. 

Fig. 2 also shows the branching ratios for different types of outcome. These are labelled 
"I" for ionization (i.e., unbinding the MS binary without triple formation), "E" for a "non- 
merging" exchange (triple formation with no star approaching within R e of the BMP), with 
either of the two MS stars ejected, and "M" for a "merger" (at least one star approaching 
within R e of the BMP). For two 0.5 M & MS stars, type M dominates when oms ^ 8 AU, and 
type I dominates for oms ^ 25AU. The formation of a triple system like PSR B1620— 26 , 
corresponding to an outcome of type E, is dominant in the interval 8 AU <^ clms ^ 25 AU, 
with the branching ratio reaching a maximum of about 50% for oms ~ 12 AU. 

In Figs. 3 and 4 we show the probability distributions of a 2 and e 2 , the semi-major axis 
and eccentricity of the outer orbit in the triple system for several typical cases. Clearly, our 
mechanism is capable of forming triples with the desired characteristics: a circular inner 
orbit (that remains unperturbed after the interaction) and a very eccentric outer orbit with 
a 2 ~ 10 2 AU. We find that the average {a 2 )/aMS ~ 3 — 5, but with a wide distribution 
extending over more than an order of magnitude. The eccentricity distribution peaks at 
some value e 2 ~ 0.6 for tight incoming binaries with oms ^ 5AU. For very wide incoming 
binaries with oms > 20 AU, the distribution has a sharp peak at e 2 = 1. For intermediate 
values of omSj the distribution of e 2 is close to thermal (and this is independent of the 
assumed distribution of eccentricities for the incoming MS binaries; see Fig. 4). 

4. Timescales and Triple Formation Rate 

Clearly, our scenario requires the presence in M4 of a significant population of wide 
primordial binaries with separations in the range 5 — 50 AU. For comparison, the maximum 
separation for a hard binary is a so f t = Gm/2a 2 = 10 A\]m ^a^ 2 , where m .5 is the average 
stellar mass in units of 0.5M 9 and a 5 is the 1-D velocity dispersion in units of 5kms _1 . The 
most recent measurements of the velocity dispersion in M4 (Peterson, Rees, & Cudworth 
1994) give a projected velocity dispersion a c % 3.5 km s -1 in the core, decreasing to ah ~ 
3kms _1 at the half-mass radius. Thus primordial binaries with separations as large as 
25 AU may in fact be (marginally) hard during most of their secular evolution outside the 
core. Even somewhat softer primordial binaries are not expected to be disrupted until after 
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they have drifted into the cluster core. This is because the characteristic time for a binary to 
drift into the core is the two-body relaxation time, which is typically larger than the binary 
disruption time by a factor of order the Coulomb logarithm In A % 10 for a marginally soft 
binary. Using eq. (5.12) of Hut & Bahcall (1983) for the binary disruption rate, we find that 
the ratio of binary disruption time to relaxation time is given (to within % 20%) by 6 



~ 100 AU may be able to drift into the cluster core before they are disrupted. 

The rate at which these wide primordial binaries drift into the cluster core today is 
~ Nb/td, where Nb is the number of wide primordial binaries still present in the outer 
region of the cluster and td ~ 10 lo yr (cf. Hut, McMillan, & Romani 1992a). We can write 
Nb = fbN Q where fb ~ 0.1 is the binary fraction and N Q ~ 10 5 . Soon after entering the 
core, the binary is disrupted in an interaction with another star or hard binary. A fraction 
/3t <C 1 of these interactions is with a (detectable) BMP and leads to the formation of a 
system like PSR B1620— 26 . Crudely we can say that /3t <^ Nbmp /N c , where ATbmp is 
the number of detectable BMPs in the core (clearly ATbmp = 1 in M4) and N c is the total 
number of objects in the core (N c ^ 10 3 for M4). We can write the average number Nt of 
triples containing a detectable BMP at any given time as Nt ~ {faNb/td^T, where is the 
average lifetime of such triples in the core 7 . Using again eq. (5.12) of Hut & Bahcall (1983) 
for a density p = 10 4 M & pc~ 3 , average stellar mass m = 1 M@, and 1-D velocity dispersion 
a = Skins" 1 , we get i T « 2 x 10 8 yr J F(£'), giving t T ~ 10 7 — 10 9 yr for £ ~ 0.1—1, with 
£ % 1.5 (a2/10 AU) -1 . Keeping all factors we obtain the following estimate: 



6 Eq. (1) takes into account only the disruption of the binary in a single strong interaction. The cumulative 
effect of distant encounters is not taken into account and may slightly reduce the value of a cr a. For very 
soft binaries (£ <C 1), this reduction would be by a factor In A 10, so that a cr ;t ~ «soft in that case (cf. 
Binney & Tremaine 1987, §8.4). However, for marginally soft binaries with £ <] 1, the results of Hut (1983) 
indicate that the reduction is by a factor of 2 at most. 

7 Note that the typical recoil velocity of the triple is negligible here and so the triple is not expected to 
spend any time far outside the core after it forms. The position of PSR B1620— 26 is projected inside, but 
most likely just outside, the core of M4. 
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Given the large uncertainties, this estimate seems perfectly compatible with the detection 
of just one triple system in M4. In addition, we note that there are now ~ 10 BMPs known 
in all globular clusters, but only one (PSR B1620— 26 ) appears to be in a triple. 
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work was begun. F. A. R. is supported by a Hubble Fellowship, funded by NASA through 
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by NASA Grant NAGW-2559 and NSF Grant AST-9308005. 
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Fig. 1. — Average eccentricity induced in the orbit of the binary millisecond pulsar (BMP) 
after an interaction with a passing star of mass 0.5 M & with a distance of closest approach r p 
to the BMP's center of mass (solid square dots). The average is over all relative inclinations, 
orientations, and phases (assumed random). The relative velocity at infinity has a fixed value 
Vgo ~ 8kms _1 . The error bars show the spread between the 10th and 90th percentiles in the 
distributions. A typical distribution, corresponding to r v ja\ = 3, is shown in the insert. For 
comparison, we also show the average eccentricity induced by a passing star of mass O.2M 
(open square dots) and 0.8 M & (open triangles). The present value of the BMP's eccentricity, 
e\ = 0.025, is shown by the dashed lines. Interactions with r p < R e % 3«i would produce an 
eccentricity larger than observed today. 



Fig. 2. — (a) Total cross section for the interaction between the BMP and a MS binary with 
semi-major axis oms containing two 0.5M 9 stars. The "effective radius" of the BMP is 
R e = 3ai, where a\ is the BMP's semi-major axis, (b) Branching ratios for various types of 
outcome: 1= ionization; E= exchange; M= merger (see text for details). The thinner lines 
show how the results vary when some of the parameters are changed: increasing R e to 4«i 
(dotted lines), increasing one of the MS masses to 0.8M 9 (short-dashed lines), or decreasing 
it to 0.2M 9 (long-dashed lines). The formation of a triple system like PSR B1620— 26 , 
corresponding here to an outcome of type E, dominates the cross section for oms ~ 10 AU. 



Fig. 3. — Distribution of the semi-major axis a 2 of the triple's outer orbit (above) and 
distribution in the (a 2) e 2 ) plane (below). The size of each square dot is proportional to the 
differential cross-section for producing a triple system with specific orbital parameters. Here 
the incoming MS binary contained two O.5M stars with semi-major axis oms — 13 AU. 



Fig. 4. — Distribution of the eccentricity e 2 of the triple's outer orbit. The incoming MS bi- 
nary contained two 0.5 M & stars. Three values of the semi-major axis are shown, oms — 3 AU 
(a), 13 AU (b), and 50 AU (c). For reference, the thin straight line shows a thermal eccen- 
tricity distribution p(e) = 2e. The eccentricity of the incoming MS binaries was assumed 
to follow a thermal distribution in all cases, except (b) where we also show the result corre- 
sponding to circular incoming binaries (dotted line). 



